Evidence for the weak steric hindrance scenario in the 
supercooled-state reorientational dynamics 
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We use molecular-dynamics computer simulations to study the translational and reorientational 
dynamics of a glass-forming liquid of dumbbells. For sufficiently elongated molecules the standard 
strong steric hindrance scenario for the rotational dynamics is found. However, for small elongations 
we find a different scenario - the weak steric hindrance scenario - caused by a new type of glass 
transition in which the orientational dynamics of the molecule's axis undergoes a dynamical tran- 
sition with a continuous increase of the non-ergodicity parameter. These results are in agreement 
with the theoretical predictions by the mode-coupling theory for the glass transition. 
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Systems of particles interacting only via hard core po- 
tentials show, in spite of their apparent simplicity, a vari- 
ety of phases: fluids, crystals, and glasses. Several inves- 
tigations have focused on the transition between these 
phases upon a change of the packing fraction tp. Re- 
cent experimental and computational studies on systems 
with elongated particles, like ellipsoids [l| and sphero- 
cylinders |2|, have revealed interesting features concern- 
ing their random (amorphous) jammed packing. E.g. it 
has been found that for the latter systems the maximum 
attainable packing fraction </? max is a non-monotonous 
function of the elongation parameter C Q which char- 
acterizes the shape of the particles Starting, for 
£ = 0, at <£> max « 0.64, the value for random close pack- 
ing of hard spheres, <y9 max reaches a maximum > 0.70 
at ( ss 0.5 and then decreases again for larger £. For 
hard core particles the quantity Aip = tp max — <p is the 
relevant parameter that governs the relaxation dynam- 
ics of the system. Hence it can be expected that the 
non-monotonous behavior of (/J max has an important im- 
plication on the tp- and ^-dependence of the dynamics, 
i.e. on the diffusion constant or the cv-relaxation time. 
For instance, it implies that, if the shape of the particle 
is changed at fixed <p, the dynamics can accelerate or be- 
come slower. Since for small Aip the dynamics depends 
very strongly on this difference, the relaxation times of 
the system can change by orders of magnitude upon a 
change of £ of just a few percent, which thus can have 
important consequences for the mechanical properties of 
the jammed packing or the crystallization dynamics. 

From the theoretical side it is possible to under- 
stand the observed non-monotonous behavior of y> m ax(C) 
at least qualitatively within the framework of mode- 
coupling theory (MCT). For the case of ellipsoids or 
dumbbells the theory predicts that the critical packing 



fraction ip c (C), at which the system undergoes a dynam- 
ical arrest, shows a maximum at intermediate values of 
C (21 0] ■ Although MCT predictions concern the line 
ific(C) an d not ^max(C)i it can be expected that these 
two curves track each other closely and that therefore 
the theory does indeed provide an explanation for the 
experimental findings. 

MCT also predicts that features of the reorientational 
dynamics strongly depend on £. It is predicted, if (, 
exceeds a certain threshold, £ > £ c (£ c = 0.345 for 
hard dumbbells |5|), that a dynamic arrest takes place 
at tpc(£) which is usually denoted as "type-B transi- 
tion" 6j. In supercooled states near <p c (()i orienta- 
tional correlation functions for angular-momentum in- 
dex £, C £ (s) (i) = (P e [e(t) ■ e(0)]), show a two-step de- 
cay with a finite plateau height (here Pi is the Legendre 
polynomial of order I and e(t) is the orientation of a 
molecule at time t.) When £ is significantly larger than 
£ c , the height of this plateau is predicted to be larger, 
the relaxation time to be longer, and the final a-decay 
to be less stretched for t = 1 than for 1=2. Such 
a scenario is called "strong steric hindrance scenario", 
and all these features are in qualitative agreement with 
the experimental results found by dielectric- loss {I = 1) 
and depolarized-light-scattering (£ — 2) spectroscopy for 
propylene carbonate |jj and glycerol [8(, and also with 
simulation results for water |i| and OTP 01 ■ 

When particles are more spherical, i.e. ( < £ c , MCT 
predicts a novel different type of dynamical transition for 
odd-£ reorientational correlators |4j, |5J, [llj . Translational 
and even-€ reorientational dynamics are predicted to un- 
dergo a dynamical arrest at a packing fraction y c (C)> 
whereas odd-^ correlators remain ergodic. Only at a 
higher packing fraction ^a(C) > Vc(C) the latter are 
predicted to undergo a dynamical arrest, i.e. the ori- 
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entation of the molecule's axis (£ = 1) freezes into a 
random direction. Most remarkably it is predicted that 
this second transition is of "type- A", i.e. the height 
of the plateau in the odd-£ reorientational correlator is 
a smooth function of ip, zero for ip < <£>a(C) and m ~ 
creases linearly for cp > ^a(C)- This is in contrast to 
the translational and even-i? rotational correlators which 
exhibit at <p c (0 a type-B transition, i.e. the height of 
the plateau already takes a finite value at (p c (()- The 
scenario for the reorientational dynamics of small elon- 
gations which is caused by the nearby type-A singularity 
has been termed "weak steric hindrance scenario" since, 
in the range <p c {Q < <P < </?a(C)> the odd-£ correlators 
relax even if the translational and even-£ reorientational 
dynamics are frozen [Tlj . In this case, the a- relaxation 
dynamics is predicted to be strongly modified for the 
I = 1 reorientational dynamics, whereas no such mod- 
ification is predicted for I — 2. So far the predicted 
existence of this type-A transition has not been investi- 
gated neither by experiments nor computer simulations. 
The goal of the present Letter is to explore whether this 
novel transition really exists and whether its features are 
in agreement with the theoretical predictions. 

Since one component systems are very prone to crystal- 
lization we have simulated a binary system of dumbbells. 
For the sake of computational efficiency these dumbbells 
are not hard core particles but instead interact with a 
Lennard-Jones (LJ) potential. Thus the system we con- 
sider is a binary mixture of rigid, symmetric dumbbell 
molecules, to be denoted by AA and BB dumbbells. 
Each molecule consists of two identical fused LJ parti- 
cles of type A or B having the same mass m, and their 
bond lengths are denoted by Iaa and Ibb- The inter- 
action between two molecules is given by the sum of 
the LJ interactions between the four constituent sites, 
V a p{r) = 4e Q/3 {(cr Q/3 /r) 12 - (cr Q/3 /r) 6 } + C Q/3 + D a(3 r, 
where a, (3 G {A, B}, with the LJ parameters caa = LO, 
oaa = 1.0, cab = 1-5, oab = 0-8, e B B = 0.5, and 
®bb = 0.88, taken from Ref. [12j. C a p and D af 3 are 
chosen so that V a p(r) and V^g(r) are zero at the cutoff 
r cut = 2.5 cr Q/ 3, and the LJ parameters are simultaneously 
renormalized to ensure that the minimum of the poten- 
tial is — e Q /3 at r = 2 1 / 6 cr ct/ 3. Bond lengths are specified 
by a single elongation parameter £, defined such that 
C = Iaa/paa = Ibb/&bb- The number of AA and BB 
dumbbells is 800 and 200, respectively, giving the num- 
ber densities paa = 800/L 3 and pbb = 200/L 3 , where 
L denotes the length of the cubic box used in the simu- 
lation. The "packing fraction" is then defined by <£>tot = 
PAA + P>bb, where p aa = (ir/6)p aa al a (l + §£ - ^C 3 ) 
for < C < 1. Fixing p t ot specifies thus L for each 
£. In this study we have used the value <p to t — 0.708. 
We note that, when specialized to ( = 0, our system re- 
duces to the spherical system studied in Ref. |l2j, apart 
from a factor of 2 in temperature T and mass. In the 
following, all quantities are expressed in reduced units 
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FIG. 1: Symbols, connected by lines, show the center of mass 
iso-diffusivity curves for the AA dumbbells. The diffusion 
constants D axe D = 1.0 x 10~ 2 , 3.2 x 10~ 3 , 1.0 x 10~ 3 , 
3.2 x 10~ 4 , 1.0 x 10 -4 , and 3.2 x 10~ 5 from top to bottom. 



with the unit of length <jaa, the unit of energy caa 
(setting ks — 1), and the unit of time (mcr AA / '(.aa) 1 ^ ' ■ 
Standard molecular-dynamics simulations have been per- 
formed like in Ref. [12j. The longest runs were 10 8 time 
steps and in order to improve the statistics of the results 
we have averaged over 16 independent runs. More details 
will be presented elsewhere 

It is necessary to make a comment on the notion of 
"packing fraction" . Since LJ interactions are finite, the 
LJ dumbbells can penetrate each other making thus the 
concept of packing fraction somewhat ill-defined. How- 
ever, it is reasonable to assume that, at fixed T, this 
penetration is only a very weak function of the elonga- 
tion (. Therefore it does make sense to study the relax- 
ation dynamics of the system at given T and t^tot as a 
function of ( and hence to check for the presence of a 
non-monotonous ^-dependence of the diffusion constant 
and the predicted type-A transition. 

To characterize the C-dependence of the translational 
dynamics we have calculated the mean-squared displace- 
ment of the center of mass of the molecules, from which 
the diffusion constant D is determined. In Fig. ^ we 
show iso-diffusivity loci, at fixed y to t, as a function of £. 
Since, for D — * 0, the iso-diffusivity lines approach the 
glass transition locus, they provide a numerical estimate 
of the shape and location of the glass line, which can be 
compared with theoretical predictions. 

From this figure we recognize that at fixed T the diffu- 
sion constant increases with increasing £, reaches a max- 
imum at around Q = 0.5 and then decreases again. Since 
this behavior is observed for all values of D we can con- 
clude that the temperature of dynamic arrest does in- 
deed show a minimum as a function of C, in agreement 
with the experimental findings 0, 0] and the theoretical 
calculations for a system of hard dumbbells (HDS) 
In Ref. |5| it has been argued that such a ^-dependence 
can be understood from the C-dependence of the first 
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sharp diffraction peak in the static structure factor and 
in Ref. ^jjj we will demonstrate that the proposed mech- 
anism seems indeed relevant for the present system. 

We now turn our attention to the ^-dependence of the 
rotational dynamics that we will characterize in terms of 
the above defined reorientational correlators (t). The 
time dependence of these correlators is shown in Fig.|2|for 
the case I = 1, ...,4. Different panels correspond to four 
representative values of C, and the temperatures are close 
to those given by the D = 3.2 x 10~ 5 iso-diffusivity locus. 
At this D value, the lowest curve reported in Fig. the 
system is strongly supercooled ^ij . For the sake of com- 
parison we have also included, as circles, the normalized 
coherent intermediate scattering function F(q, t) for the 
wave-vector q that corresponds to the first maximum in 
the static structure factor [lflj . 

Most of these correlators show a two-step relaxation: 
They decay on intermediate times to a plateau, a motion 
that corresponds to the relaxation of the molecule inside 
the cage formed by its neighbors, followed by a final a- 
decay in which the molecules leave their cage (or reorient 
in the case of Cf ] (t)). In the following we will denote the 
plateau height by and estimate its value via appropri- 
ate fits of Cf } (t) with von Schweidler's law The 
a-relaxation time ti will be defined as C^{ti) = ff/10. 

For sufficiently large elongation (( = 0.8, Fig. 0a,), 
steric hindrance for reorientations is substantial, and thus 
cf\t) exhibit the dy namics of the strong steric hin- 
drance scenario 0, |n| : The plateau height fl and the 
a-relaxation time ti decrease monotonously with increas- 
ing £. In particular, ff > /f and n > r 2 . These features 
are in qualitative agreement withprevious findings for 
various glass-forming systems 0, la III El ■ Upon a de- 
crease of the elongation, the steric hindrance effects for 
reorientations weaken, leading to smaller compared to 
those for larger elongation. The relations ff > and 
T\ > Ti hold until £ « 0.6. For intermediate elongations, 
e.g. £ = 0.5 shown in Fig. |2h>, the condition /f > /| still 
holds, but now t\ < t 2 . Thus for intermediate elonga- 
tions one expects that the strength of the a-peak in di- 
electric spectroscopy - a technique sensitive to the 1=1 
decorrelation - is more pronounced than the one in dy- 
namic light scattering (£ = 2), but that the location of 
the former is at higher frequencies than the latter. 

If the elongation is decreased even further, Figs.0: and 
Eli, the a-relaxation process for £ — 1 basically vanishes, 
although the one for I = 2 is still clearly visible. In this 
case fi = 0. A closer inspection of the curves shows 
that this result holds for all odd values of I. Thus, for 
small C, a dielectric spectroscopy experiment would not 
show anymore an a-peak, although such a peak would 
still be visible in a dynamic light scattering experiment or 
an inelastic neutron scattering experiment (which probes 
F(q, t)). Results reported in Fig.|21provide evidence that, 
for nearly the same value of D, the relaxation dynamics of 
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FIG. 2: The coherent intermediate scattering function (cir- 
cles) and the reorientational correlators (t) for the ^4^4 
dumbbells in supercooled states for I = 1 (solid line), 2 
(dashed line), 3 (dotted line), and 4 (dash-dotted line) for 
(a) C = 0.8 and T = 1.06, (b) £ = 0.5 and T = 0.75, (c) 
C = 0.3 and T = 1.13, and (d) ( = 0.2 and T = 1.85. 



the molecules depends strongly on their shape: For elon- 
gated molecules all correlators decay on a similar time 
scale, whereas for more spherical molecules the transla- 
tional and even-^ reorientational dynamics relax much 
slower than the reorientational ones for odd values of £. 

The observed peculiar behavior of the odd-£ reorienta- 
tional correlators at small values of C can be interpreted 
as being caused by the nearby type-A transition, pre- 
dicted by MCT for hard dumbbells Evidence for the 
existence of this transition can be obtained by investi- 
gating the ^-dependence of the height of the plateaus in 
the orientational correlation functions, i.e. /|(C)- MCT 
predicts that, moving along the ideal glass line, for even 
I this height decreases smoothly with decreasing £ and 
becomes zero at C, = 0, reflecting the weakened steric 
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FIG. 3: Plateau heights (circles) and /| (filled squares) 
for AA-dumbbell's correlators C^\t) for £ = 1 and 2, respec- 
tively. The dotted line represents /f for the HDS taken from 
Ref. |fj with the arrow indicating the critical elongation C™ s . 

hindrance effects. In contrast to this, the ^-dependence 
of fgiC) for odd I is predicted to show a smooth depen- 
dence on £ down to the critical value £ c discussed above, 
to vanish at £ = £ c and then to stay zero for £ < Q c . 

These predictions are indeed compatible with the re- 
sults from our simulations as demonstrated in Fig. |3| 
where we plot the ^-dependence of fl and /f . Starting 
from large £, the plateau height /| decreases smoothly 
toward the limit /f = for spherical system (£ = 0). 
Instead, the height /f decreases smoothly with decreas- 
ing £ at large and intermediate £, but it shows, around 
Q = 0.4, a sudden drop and becomes zero within the sta- 
tistical accuracy around £ = 0.3. This behavior suggests 
the existence of a critical elongation ( c which is located 
between £ = 0.3 and 0.4. For the sake of comparison we 
have added to Fig. also the theoretical curve /f for the 
HDS (dotted line) as well as the location of its critical 
elongation C? DS |5j. Note that shows nearly the 

same ^-dependence as the theoretical quantity for HDS, 
which corroborates the existence of £ c between £ = 0.3 
and 0.4, and hence the presence of the type-A transition 
in molecular systems of small elongations. 

Since the appearance of the plateau regime for the 
even-£ correlators is due to the cage effect, the results 
from the simulation for £ = 0.2 and 0.3 (Figs. [5J: and 
Eli) imply that (i) the relaxation of the odd-£ correla- 
tors occurs via large-angle flips between the energetically 
identical orientations e and — e, which do not alter the 
even-£ correlators, and (ii) these flips occur before the 
molecule has left the cage. Such a behavior for the reori- 
entational dynamics is in striking contrast to the strong 
steric hindrance scenario for large elongations for which 
the reorientation occurs continuously and on the time 
scale of the molecule's escape from the cage. 

Even for elongations exceeding £ c , precursor effects of 
the type-A transition do influence seriously the results 
for the odd-£ reorientational dynamics, the most notice- 



able one being the speeding up of their cv-processes . 
For instance, one concludes from Fig. 03 that the steric 
hindrance effect is still considerable for ( = 0.5 since 
fi remains close to one. However, the canonical order 
Ti > T2 expected for the strong steric hindrance scenario 
is violated for ( — 0.5, as can be seen from Fig. 12b, which 
is evidence for the nearby type-A transition. 

Although strictly speaking the type-A transition dis- 
cussed here exists only for symmetric molecules in which 
two ends are identical, most of the mentioned features 
can be expected to be observable, in the form of pre- 
cursor effects, also for molecules in which this symmetry 
is slightly broken. E.g. Kammerer et al. have stud- 
ied a system with slightly asymmetric dumbbells with 
C = 0.5 0| . These authors found t\ < t 2 in spite of 
the large plateau /f , in full agreement with the present 
results and the theoretical expectation, indicating that 
the scenario caused by a nearby type-A transition can 
indeed be found in a class of molecular systems which 
are weakly elongated. 
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